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ABSTRACT 

A multilayer model of the earth is used in finite element calculations of time dependent de- 
formation and stress following an earthquake on a strike-slip fault. The model involves shear 
properties of an elastic upper lithosphere, a standard viscoelastic linear solid lower lithosphere, a 
Maxwell viscoeiastic asthenosphere and an elastic mesosphere. Elastic dilatational properties are 
assumed throughout. Time dependent displacements, strains, and stresses are computed both at 
the surface of the earth and at depth. The analysis includes both systematic variations of fault 
and layer depths and comparisons with simpler elastic lithosphere over viscoelastic asthenosphere 
calculations. For conditions which may be appropriate for the earth both the creep of the lower 
lithosphere and asthenosphere can contribute to the postseismic deformation. The magnitude of 
the deformation is enhanced by a short distance between the bottom of the fault (slip zone) and 
the top of the creep layer but is less sensitive to the thickness of the creeping layer. Further- 
more postseismic restressing is increased as the lower lithosphere becomes more viscoelastic, but 
the tendency for the width of the restressed zone to grow with time is retarded. 
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A MULTILAYER MODEL OF TIME DEPENDENT 
DEFORMA VION EOLLOWINC. AN EARTHQUAKE ON A 
STRIKE-SLIP FAULT 


INTRODUCTION 

There is r ^mwing btxly of cvidena* tiuil ix)stseisnuc givuiul deformation can be of geoiletic 
and tectonic significance following a major eartlu|uakc. In a in aiber of cases there have been 
direct observations of (x^stseismic vertical and/or horizontal deformations (e.g. Ehatchcr, 1^75. 
Brown, et al.. ’.'>77; Thatcher, et al„ I '>80, Dunbar, et al„ 1^80). There are, m fact, a number 
of mechanisms that have the {x>tential to create significant ixjstseismic deformations and identify- 
ing the operative mechanisms in specific cases may prove more difficult than detecting the 
motion. Some of the (Xitential contributors to postseismic deformation that have been suggestal 
include |x>st«eisinic fault slip, visanMastic relaxation of subsurface layers of the earth (Nur and 
M.ivko, l'>74; Melosh, l'>70; Rundle and Jackson, I '>7V , 1 hatcher and Ruiulle, l‘>7*); (’ohen, 
l')SlaK fault creep at depth (Thatcher. I‘>74j, How of low visaisity magma regions iWalir aiul 
Wyss. I')S0), and other forms of anelastic relaxation of surface and subsurface material (('ohen, 

I '>80 a, b; Yang and I'okso/, I '>8 IT I'urthernHue the pattern of deformation can be intluenced 
by horizontal and vertical variations in elastic (Mahrer and Nur. I'>7*)J and anelastic material ci'n- 
stants, coseisinic slip orientati* ns (Mansinha and Sinylic, '>7n, and fault length, width, depth, etc. 

Ihe mechanism of crustal deformation of interest in this study is visaiclastic How of sub- 
ciustal rocks. It is well e.stablished by ex|K'nmeiits in rock mechanic's that anelastic creep is a 
function of nx'k type, temperature, stres;;. pressuiv and |x*rhaps other parameters (('alter, I'^^tr). 
There is little doubt that there are conditions in the earth that can lead to some form of an- 
el.istic creep (.Stocker and Ashby. 1973; (ieot/e and Evans. 1 979). therefore it is not unreasmiable 
that viscoelastic res|xnises have been invoked in studies of anivection. |xistglacial rebound, earth 
tidal detoi Illation, seismic wave alteiuiatiun, and more receiitl> |x>stM'isjiiic relxniiul, Elie foiin 
of Ihe viscoelastic cxuistitutive law is somewhat less certain. I'heie are a nrmber ol laborator> 
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experiments for winch non-tincer How bws. of the form < ■ A iH' f(P/T) where « is the iteviutoric 
strain rate, o is the deviatoric stress, f is a function of pressure, P, and temiieratuiT, T, and n — 3 have 
been applied with success tWcerlman and Weertman. h>75l. It would m>t be surprising! if a simi- 
lar law is o|w?rative within some regions of the earth; however, the lack of kimwledge of the 
earth's interior state, the existence of non-steady state stressing, and the existence of other geo- 
physical complexities make the form of the rhei>logical law or laws uncertain. In th>* present 
analysis ! intend to use simpler linear ^isciu'lastic rheologies in the development of nutdels of 
(vstseisiuic deforinatmn. The sclectmn of a linear iiuhIoI is a simplification, bul one which 
allows the separation of the a>s«*ismi> _.id pi>stsoisi,uc delorm,itk>n and stress Helds fr»>in pre- 
setsmic ones. Although the nnHlel is ci'nci'ptiullv clear, the mathematical ,»nd ivinputational 
detail is sufi'iciently ivmplex that a numencal. tinieAlet>endent finite ef-mient soIiiIkmi i»f the 
governing eituations is justified. 

MODI L DKSC RIPTION 

The nuHlel chosen for this study is shown in I'igure I. the rationale for the model has been 
discus.si'd in I'ohen |l*>SlaV It is a multilayrr nunlel with a strike-slip fault onibeddetl in the 
upper layer. The vertiad strike-slip fault extends from an uptvr depth which may be zero (the 
surface ol the earth), to a depth»/ . The upper layer, the uptvr lithospherv, extends fix>m the 
surface to a depth, II, (generally M, i?* 'ybut there may be exceptions). I his layer is chai.icter- 
i/.ed by the elastic parameter ncidity,p|. Below this layer to a depth llj. liC' the k'wer lithosphere 
represented by a standaid viso>elastic si>IhI with rigulitics and vtscosily. »j> as sliowii 

in the figure. Below the lithosphere, to ,i depth. lies the aslhenospheie reprevnled by a 
Maxwell viscrvlast ic lluid with rigidity, ity and viscosity, tj. Below the asthenosphere is an 
clastic mesosphere with righlity. The aforementk'neU rhixvlogieal im»dels n*pn*\ent the earth's 
deviatoric pro|'»erties. the dilatatkmal pa>(H.Ttk's are avsiimed elastic (ImMighout with avw>ciale»l 
bulk nuxluli k,, kj. kj. am' k^. The clH*ice of a standard linear si'lid k>wer lill.osphere allows 
for .«n inituil elastic response (rigidity, p,) to an applievi constant strain, a iKiituI strevs leLxation 



(on a tunc scale dctcrinincd by and the rigidities), and permanent support of a residual strv'ss 
with a reduced rigidity, p, Similarly, the Maxwell viscoelastic t 1 uid model ol' the 

asthenosphere allows an initial elastic respimse (rigidity, p^) to an .ipplied constant strain and 
subseiiuent stress relaxation (on a time scale controlled by t/j and p^). 


The mathematical description of the model can be developeil by considering a stbslance 
that has standard linear solid deviatoric and elastic dilatational properties. This is the lower lith- 
osphere element of the problem; other layers will be derived by evaluating various limiting forms 
of this general substance. The standard linear solid consists of a parallel annbination of elastic 
and viscous elements (Kelvin element) in senes with an clastic lenient. Although only the et|ua- 
lions for the shear components of stress and strain, e.g. f j j and e will be needed in this analysis, 
the equations lor all ci>m|'onents will be derived as they will be needed for dip-slip analyses 
(Co/cn. I^SIb) and three dimensional calculations. I'or the elastic clement of the Kelvin clement 
the deviatoric constitutive relations are of the form 

20ii - O22 -O33 (2t,| - t22 -t 3 j) 

_ _ _ = I 1 ^ 


“12 “ ^12 


(lb) 


wliereO|i and f,i are. for example, the ( I . I ) aunponents of stress and strain respectively, and the 
quanlilies t2o|, - Oii - 033)/.^ aiul (2c.,, - t ,, - t are the corres|xinding deviatoric stresses 
and strains ol the normal eoniix>nenls. Tor the viscous part of the Kelvin element the con- 
st it iilive relations are 


2o , , -1)22-0 3, 

3 



(2 a) 


(1,2 = 2fj 


dt 


(2b) 


.Since llie Kelvin elements are connected in parallel the strains in equations 1 and 2 are equal and 
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the stresses add. tluis for the Kelvin element 


(3a) 


2o,, - 022 “ nj 3 « 2 ♦ r? — ^ (2«,, - ejj “ ^ J3^ 

Oi2 ’ ' '>2 


(3b) 


or formally wiih 
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12 


'12 


2{fi^ + j?D) 

.In analogy vwitli equation I, for the elastic clement in series with the Kelvin substance. 

2 o,| - O 22 - O 33 


(4a) 


(4b) 


-f,i - €22 - fi 3 = 


‘12 “ 




0|2 


(5a) 

(5b) 


The series elements have equal stresses and their stniins add, hence for the combined three ele- 
ment standard linear solkl 


'■ - [»'i Mb * 'i»] 


((ni) 


((>b) 


or rearranging 


2 o|, 03 J + 


1^.1 *’ 


(20|| - fl;j - Ojj» 


- ^22 ■ ^33' ■* (-^11 ■ ^22 ~ *33' 

n 


U 12 I 


(^.1 ♦ Ub) 


0|J = 2u^ e,, + 


^b 


f|2 


(7a) 


(7b) 


when the dots above the svm'.H.Is indicate time derivative;. 
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The next step is to impose the condition that the dilatation is elastic 


0|1 •*•022 + 0,3 * 3 k ( f ,, + 622 + C „) 


Siilving equation 8 for O 22 + o,, and inserting in equation 7a results in 
Oia+Mb) r 3 k + 4 #i 3 l 

‘Ml + o,, = I I e„ + I — (622 + ^33! 


l^b ■ 

3k + 4 

Aig+#ii, 

I? 
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^11 ■*• 


„ '"a Alb 

3k - 2 ( Ma + Mh > 

Ala + Alb (f 22 •*■ ^ 33 > 

ff 
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«; • . Al, Alb . 

Writing - and r = + /ij,) the general txiuations are 

Ala •‘■Alb 


0,1 /3k + 4fi\ /3k - 2^\ /3k + 4ju'\ ej, 

_ = _ . __ . 


3k - 2n’\ / 6,2 + f33 


O 22 /3k + 4^\ /3k - 2ti\ /3k + 4/a’\ e-,-, 

/3k - 2fi'\/^^^ + 6 „\ 


03 , /3k + 4^\ /3k - 2#i\ /3k + A^l'\ e,, 

/3k - /t,, + 622 \ 


0,2 + = ^/i 6,2 + 6,2 


(ll)d) 


= ./i 6 |, + e„ 


= .^ 62 , + 6,3 
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Equations 10 form a system of equations with standard linear solid de' iatoric and clastic diluta- 
tional properties. A system with Maxw«..l dcviatoric properties can he derived by letting Mb 0. 
This implies ft' 0 and, fur example 


/3k ♦ 4 m\ /3k - 2it\ I 

♦ '33' - 7 


0| I" k (e 1 1 e->7 e ^i)l 


II 


33' 


(lla) 


d,j + — • 2Mif|: (lib) 

r 

A purely elastic substance can he derived from equation 1 1 by letting r tlicn 

/3k + 4 m\ /3k - 2 m\ 

°ii “ ( — ^ — 1 ^ ( — ~ — j (‘22 ^ 33 ! 

In the finite element analy~is I have used a two dimensional grid in wh:,n the directions are 
depth and distance from the fault. The displ icements (for a strike slip fault) arc orthogonal to 
this plane in the direction of the 3 axis). FIk' technique for solving etpiations 10 in the finite 
element analysis is presented in Appendix A. 


Having derived the rheological ec|uations. and implemented the finite element technique for 
solving them, the next step for computing the quasi-static displacement, strains, ami stresses is to 
define the input numerical parameters of the model. The parameters that affect the model calcu- 
lations arc the fault slip, the depths of top and bottom lH3undari(*2> of the f alt slip rone, the 
layer boundaries H|, Hj. and Hj and the clastic and viscous p;irameters of the layers. For sim- 
plicity I have cho.scn to set equal all the rigidities with Mi ■* M 21 * M^h * * fa * X 10" 

dyne/cm^ and bulk nK>duli with k| * kj • kj » k 4 » 8.33 X 10" dyne/ciiW. These values cor- 
respond to Young's moduli of 1.2.S X lO'^dyne/cm^, and Pi'isaon’s ratio ol 0.25. The viscxvsiiy 
of the lower lithosphere is chosen to be I X lO^l poise and that of the presumably wanner astlu no- 
sphere 5 X 10*’ |ioisc. The .sssociated time c'onsiants for stress relaxation are * 
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^2b) “ i X 10’ sec — 32 yrs., and T3 ■ 1 XIO* sec » 3 yrs. Frequcnily in succoedint! 

paragraphs I will examine results at five particular times: 

tj, = O’*" sec 

tg = I -10* sec *= = 0,1 T2 ^3 yrs 

tj, = 5 ■ 10* sec * Stj * 0.5 f2 16 yrs 

tj. * 1 • 10’ sec * I0 t3 = T2 *■32 yn 

tj = 5 • 10’ sec = 50t. * 5 t 2 ** 159 yrs 

The slip function is set equal to I meter from the top of the fault (the surface of the earth in 
the present calculations) to a depth D; then by way of continuity of displacement within the 
finite elements displacement decreases linearly toward 0 at the node below D (generally 10 km 
lower). Ihe slip enters the calculations as a scale constant: thus, the repoied values of displace- 
ments. strains, and stres.ses arc normalized per one meter slip. The bottom of the computational 
grid is arbitrarily s«*lcdcd at XOOkm. the ends of the grid at ±4000kin sufficiently far from Ihe 
source region that the boundary conditions (free or fixed) do not affect the results. The remain- 
ing parameters of the model arc D, and li|, the depth of the boundary between the upper and 
lower lithosphere, H2. the depth of the lower lithosphcre-asthenosphc.e Ixiundary and H3. the 
depth of the a.sthenosphenr-mesosphere boundary. I designate model calculations by the notation 
(D), ll|, M2. II3 where the numerical values are given in kilometers. 

Rism.s 

Ihe starting point for presenting the results of model calculations will be model (20). .10. 

75. 4()() for which Ihe cxrseismic slip p’ane terminates at the top of the lower lithosphere at a 
depth of thirty kilometers. The calculated displacemc.its are shown as a function of distance 
ln>m the fault and time followh'g the earthquake in I'igure 2. Appieciahle ptisiscismic displacc*- 
inents extend over a broad /one extending to scver.il huiulicd kilome'ers fiom the fault. The 
jv.ik displacement is about 7cm lO time tj, !(> yrs. and about 15cm at time tj 159 yrs. 

There are considerable similarities in the displacement versus distance curves at various times. 
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howt*ver, the location of the peak displacement does vary somewhat. This peak point, which is 
a tram \ ion pi>int from positive to t entative postscismic strain, C 13 , is at X — 150km at time t^ ~ 

3 yrs. and at X — 200km at t^ ~ IS9 yrs. Tlic details of the time dependent motion of the 
peak will be discussed below. 

It B instructive as to the development ol the postseisinic deformation to exami e the sub- 
surface displacement and stress fiekts. Figure 3 shows the coseismic displacements at three depths: 
the surface of the earth, the bottom of the slip region, (i.e. the top of the lower lithosphere), and lop 
of the asthenospherc. The gradient of these curves gives the coseismic strain, 6 , 3 , (in all the rru- 
mcrical results presented herein engineering strains are used) or to within a multiplicative constant, 
the initial stress. O 13 , acting in u plane pirallcl to the fault. The curves show the initial stress 
drop (negative gradient) at the surface of the earth and significant stress rises ()>ositive gradients) 
near the fault below the slip region. In response to the applied stresses the viscoelastic regions 
begin to flow. The flow produces a time dependent stress relaxation as shown in Figure 4a. The 
asthenoNphere, which has a shorter relaxation lime than the lower lithosphere, accomplishes alxiu* 
9091 of its a>mplete. near fault, relaxation of the coseismic stress of 0.51 X 10* dync/cm^ within 
about 5X10* sec. ?: 16 yrs. The more sluggish lower lithosphere responds over a long time scale 
and of course maintains a r sklual permanent stn.*ss field. There is significant redressing 
of the surface that accompanies the subsurface relaxation and flow. TIk time dependent 
surface ami subsurfao* strain changes accompanying these stress changes arc shown in 
Figuie 4b. The detailed dependence of the surface stress and strain on distance from the 
fault and time is siiown in Figure S. Althougli the postscismic strains are ciiinmonly 
*inailer m magnitude than the coseismic ones, they are more broadly distmuted. For 
example at time, 1 ^, the postscismic stress is reduced by om'-lialf fix in its value over a 
distance of about 35km; by contrast the aiscismic strain drops by one-hall over 20km. The 
broad rone of postsetsmic deformation is. of course, a reflectKin of the deep |m>siIioii of the 
tViwing viscoelastic layers. 
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Another way of studying the spatiat and temporal dependence of the stress (or deformation) 
fields is to plot contouis of constant postscismic stress on a distnnee-tinie yrid as ::hown in l-igiire 
6a. Th's sketch illustrates that the region of positive postseisinic stress advances from within 
^135 km at t ■ I X 10* sec to -*-189 km at t ■ .^X 10’ sec. Other contours of constant strc'ss 
show similar progressions away from the fault with time. Figure 6b shows the corresponding 
contours of constant stress when the lithosphere is modeled as an clastic substance rather than as 
clastic layer over a viscoelastic one. Particularly noteworthy in Figure 6b is the greater extent 
of the positivi* |X)stsei»mic stress zone, noticeably at long times. The damped advancement of 
the postscismic stress field in the case of the layered viscoelastic profile can be understood by 
considering the depths and response times of the layers. First as the asthenosphere relaxes it 
inducc*s surface displacements which are zero at the fault, rise to a |X'ak. then decrease further 
away. As relaxation progresses the peak moves further from the fault implying an advancement 
of the Z' ne of postscismic restressing. As the lower lithosphere begins to relax, however, it too 
produces si'.rtace displacements, ones that will be relatively near the fault. Thus the advance of 
the displace.nent peak away from the fault with time will be retarded and the zone of positive 
restressing iianowed, but the magnitudes of the near fault stress recovery will be increased. 1 his 
|H)int IS emphasi/cd in Figure 7 win h shows the rate of horizontal displacement versus distance 
and time. A comparison of Figure 7a for the viscoelastic lower lithosphere model and Figure 7b 
for the elastic lower lithosphere models sliows that in the former case and at long times there is 
a hroadei zone of significant displacement rates. 

rile eoseismic and postscismic displacements vary not only with distance uom the fauli. but 
also with depth from the surtace. This gives use to the stress and strain components, Oj). and 
t ris|H-ctively. Whereas d|j acts on a plane parallel to the fault, a acts on a plane parallel *o 
the earth's surtace. As Figure H illustrates is generally less than o^y at most points near the 
surtace but bcct»mes of greater relative significana* at depth. The cxirresponding depth, distance, 
and iiiiK dependent strains are shown in Figure 9. 


Ai iUuft.'^teti carhei H it iratruiiive to cnmparr the pment culculttmns of surface def o rnw 
tiont »*fth Ihoic which ignore vncocbsticity in t'le k<wer lithosphere hy citnsidcrinit an ctaalic 
lithosphere over a viscoebstic asthenospherc. Figure 10 shows a comparison of the calcubted 
dhph cements s'srsus dntanu' and rime. The cuirrs are sitniUir in shape in the two modeb but 
differ significantly in magnitude, most ruitii-cably at long time when relaxation of the upper 
lithosphere has had sufficient opportunity to occur; at t ■ t^. the maximum difference ui the 
curves is about 6 cm at X ■ 75 km. Of pertups greater significance is the diftcrencc in strains. 

Figure 1 1 shows the ratio of the nusleled strains versus distance and time. Near the fault the 
effect of creep in the lower lithosphae b to produce a three-fold increase in the postscismic 
restraining over a time of I years. 

Another aspect of the calcubtiun that henrs investigation b the effect uf' > hangcs in the layer 
depths on si*rfacc displacements. Figure 12a, for example, shows :hc effect of <. hanging the depth 
of the lower lithosphere- jsthenosphcrc l«oundary. Hj. after a time tj, w 16 yrs. Of particubr note 
is the fact that the position of the peak of the postveismic dispbcetnenl pattern is fairly sensitive 
to the depth of the top of the asthenosphere. The ck>ser Hj is to the bottom of the slip reg»n 
the greater the initbl stressing of the ast her- sphere and the greatn the subsequent p.ist seismic 
defonration. Fhe sensitmty of the deformaikin pattern to H 2 suggests that observations of 
postseisniK dispttcemcnts may be useful for determining the depth to the top of the asthenosphere 
(if such obtervatior.s can he sufliciently stripped of other amiaminating s^tnals fram. e.f.. plate 
motion). By contrast tlir location of the bsrltom of the asthenc'';phcre b not well resolved by 
the r.hape of the displacement curve or b> ilie amplitude (Figure 12b). Similarly the shape of the 
dirpiacement vs. distance curve at tiir>e tj is iku very wnsitiv.. to tiir depth of the top of the 
tTeepiigi lower lithosphere (Figure 12c) although, the dbpiaerments ar<* enhanced as the distance 
between titc creep rone and the slip area u reduced. Tfu Ncnsitivity of the displacement curve 
shape to H| b increased only slighlly when the displacetrH-nls at time t), are subtracted fiom 
those at tj to estimatrthc long time scale conip<>neni 01 'tisplaccment. The effect ol leducinr 
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fault lower depth on coseismic and postseismic displacements are shown in Figure 13. Both the 
coseismic and postseismic displacements are larger for the more deeply extending fault. However, 
while the coseismic strain drop is reduced, the near-fault postseismic strain recovery is increased 
by increasing D. 

CONCLUSIONS 

The viscoelastic rebound effects discussed in this paper will be significant when the depths of 
the anelastic relaxation layers are not much greater than the tM^ttom of the fault slip zone. There 
is insufficient data currently available to decide whether such a situation exists on any major strike- 
slip fault of the earth. Predicted peak postseismic strains of lO-lOOpstrains following major earth- 
quakes are, however, detectable by surveying techniques, although measurement interpretation may 
be ambiguous. The model calculations presented here have shown that the spatial and temporal 
patterns of postseismic deformation can be sensitive to the depth of the creep region and the 
viscosity structure of the earth, however, the spatial-temporal pattern is less sensitive to creep 
region thickness. It is possible that definitive models of the deformation of the earth during 
the earthquake cycle will require attention to multilayer viscoelastic effects as well as other fac- 
tors such as lateral and vertical variations in slip, rheological properties, rock type etc. 
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APPENDIX 


Mriosli Mnd Racfsky (1980) have discussed the use of the explicit algorithm of Cormemi 
(197$) and the implicit alforithm of Hughes and Taylor (1978) for finite element analyses cf 
equations of the form 


a * D(f - e^P ) 


(Al) 


where D is the matrix of elastic parameters and is a funetk^n of stress. The explicit algorithm 
cited above is also applicable when e''P is a function of strain as well as strevs. and the reader is 
referred to these papers for a description of the technique. In this appendix etfuations 10 of the 

text are recast in the form of equation A I so the algorithm cun be applied. Specifically applying 
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equation Al to the normal comtx>nents of a in equation 10 results in 
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It is aiuvoiuont to clnuiiwtc Ojj thioujih i«quation 8. I'hon 
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FIGURE CAPTIONS 


Figure I. Multilayer rheological model for postseismic rebound study. 

Figure 2. Displacement versus distance; coseismic displacement at time t^; postseismic displace- 
ments (AW = W(t) - W(t(,)) at times tj = 1 • 10* sec 3 yrs, tj, = 5 • 10* sec « 16 
yrs, tg = 1 • 10^ sec 31 yrs, tj = 5 • 10^ sec 1 59 yrs. 

Figure 3. Coseismic displacements versus distance at various depths. 

Figure 4a. Ntar-fault stress versus time at various depths. Stress is averaged over a grid element 
extending from 1-lOkm from the fault and over indicated depth ranges. 

Figure 4b. Near-fault strain versus time at various depths. 

Figure 5. Stress (strain) versus distance from fault at various depths. 

Figure 6. Contours of constant post.seismic stress (units of 10* dyne/cm^) on a distance from 
fault-time after earthquake grid. a. Model (20), 30, 75, 400 with viscoelastic lower 
lithosphere, b. Model (20), 75, 75, 400 with elastic lower lithosphere. 

Figure 7. Rate of horizontal displacement versus distance from fault at various times, a. Model 
(20), 30, 75. 400. b. Model (20), 75, 75, 400. 

Figure 8. Stress components, a 13 and 033, versus distance from fault at various times and depths. 

a. Depth = 0-1 0km, b. depth = 30-40km, c. depth = 75-1 00km, (note vertical scale 
change for c). 

Figure 9. Strain components. 613 and e23> versus distance from fault at various times and depths, 
a. Depth = O-IOkin, b. depth = 30-40km, c. depth = 75-1 00km. 

Figure 10. Comparisons of postseismic displacement patterns at various times for models with 
and without a viscoelastic lower lithosphere. 




Figure 1 1 . Ratio of near surface postseismic strain for viscoelastic lower lithosphere to near su^ 
face postseismic strain for elastic lower lithosphere as a function of distance from 
the fault and time. 

Figure 12. Postseismic surface displacement versus distance from fault for various layer interface 
depths and selected times, a. Lower lithosphere - asthenosphere boundary, H 2 , varied, 
t * tj,, b. asthenosphere - mesosphere boundary, H 3 , varied, t = tj,, c. upper litho- 
sphere - lower lithosphere boundary, Hj, varied, t = tj. 

Figure 13. Coseismic and postseismic displacements versus distance from fault for different 

width faults. Solid lines model (20), 30, 75, 400, dashed lines model (10), 30, 7S, 
400. 
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Fkure 2 DispUocment versus distance; ooseismic displacement at time Iq; 

(AW . W(.) - W.„)) .. I. J ^ lOj K. • 3 

tfc = 5 • 10 * sec • 16 yrs. tc = I • lO’ sec ^ 31 yrs. ta - 5 W sec 

159 yrs. 
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Figure 5. Strain (stress) versus distance from fault at various depths. 
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Figurr 6. Contours of comtani pottaeismic itre» (uniU of 10* dync/cm^) un a d»t:inLX' 
from fauh-fimr after earthquake fhd. a. model (20), 30. 7S. 40U with viMXJclaslic lower 
lithotpbers, b. model (20), 75, 75, 400 with elastic lower lithosphere. Ij 
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Figure 7. Rate of horizontal displacement versus 
model (20), 75, 75, 400. 
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Figure 12. Postseismic surface displacement versus distance from faull for various layer inter- 
face depths and selected times, a. lower lith*jsphere - aslhenosphere bouiulary. Hj, varied, t = t,,. 
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l-igurc 13. Coscismic and postscismic displacements versus distaiux frum fault lor 
different width faults. Solid lines model (20), 30, 7.S, 400, dashed lines model 
(10). .30. 7.S. 400. 
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